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Geodesically equivalent metrics in general relativity 

Vladimir S. Matveev* 



Abstract 



We discuss whether it is possible to reconstruct a metric by its unparameterized 
geodesies, and how to do it effectively. We explain why this problem is interesting for 
general relativity. We show how to understand whether all curves from a sufficiently 
i -^h ■ big family are umparameterized geodesies of a certain affine connection, and how 

to reconstruct algorithmically a generic 4-dimensional metric by its unparameter- 
ized geodesies. The algorithm works most effectively if the metric is Ricci-flat. We 
also prove that almost every metric does not allow nontrivial geodesic equivalence, 
£Nj ■ and construct all pairs of 4-dimensional geodesically equivalent metrics of Lorentz 

| signature. 

o i 

1 Introduction 

Let (M n ,g) be a connected Riemannian (= g is positive definite) or pseudo-Riemannian 
manifold of dimension n > 2. We say that a metric g on M n is geodesically equivalent to 
g, if every geodesic of g is a (reparametrized) geodesic of g. We say that they are affine 
equivalent, if their Levi-Civita connections coincide. 



The first examples of geodesically equivalent metrics are due to Lagrange [32] . He observed 
that the radial projection f(x,y,z) = (— |, — |, — l) takes geodesies of the half-sphere 
S 2 := {(x,y,z) G R 3 : x 2 + y 2 + z 2 = 1, z < 0} to the geodesies of the plane E 2 := 
{(x, y,z) G IR 3 : £ = —1}, since the geodesies of both metrics are intersection of the 
2-plane containing the point (0, 0, 0) with the surface. Later, Beltrami [31 H] generalized 
the example for the metrics of constant negative curvature, and for the pseudo-Riemannian 
metrics of constant curvature. In the example of Lagrange, he replaced the half sphere by 
the half of one of the hyperboloids H± := {(x,y,z) G M 3 : x 2 + y 2 — z 2 = ±1}, with 
the restriction of the Lorentz metrics dx 2 + dy 2 — dz 2 to it. Then, the geodesies of the 
metric are also intersections of the 2-planes containing the point (0, 0, 0) with the surface, 
and, therefore, the stereographic projection sends it to the straight lines of the appropriate 
plane. 



"Institute of Mathematics, FSU Jena, 07737 Jena Germany, vladimir.matveev@uni-jena.de 



1 



Though the examples of the Lagrange and Beltrami are two-dimensional, one can easily 
generalize them for every dimension (for Riemannian metrics, it was done already in [3]) 
and for every signature. 

Since the time of Hermann Weyl, geodesically equivalent metrics were actively discussed in 
the realm of geneal relativity theory. The context of general relativity poses the following 
restrictions: the dimension is 4, the metrics are pseudo- Riemannian of Lorentz signature 
(—,+,+,+) or (+,—,—,—), and sometimes the metrics satisfy additional assumptions 
such that one or both metrics are Ricci-flat (Rij = 0), or Einstein (R^ = jgij), or, more 
generally, satisfy the Einstein equation — ^g%j = Tij with 'physically interesting' stress- 
energy tensor Ty. 

Let us explain (using a slightly naive language) one of the possible motivations for this 
interest. Suppose we would like to understand the structure of the space-time in a certain 
part of the universe. We assume that this part is far enough so the we can use only 
telescopes (in particular we can not send a space ship there). We still assume that the 
telescopes can see sufficiently many objects in this part of universe. Then, if the relativistic 
effects are not negligible (that happens for example if the objects in this part of space 
time are sufficiently fast or if this region of the universe is big enough), we obtain as a 
rule the world lines of the objects as unparameterized curves. Indeed, local coordinates 
on a 4-manifold are 4 smooth functions on the manifold such that their differentials are 
linearly independent. Now, for every freely falling object in this part of the universe such 
that it can be registered by telescopes, each telescope at every moment of time gives us 
two such functions, namely the spherical coordinates 4> and 9 (latitude and longitude) of 
the direction the light reflected from the object comes to the telescope from (in a naive 
language, the telescope 'sees' the direction where the object lies), see the picture below. 
Since we have two telescopes, altogether we have 4 functions of t, (4>i(t), 6i(t), (j>2 (t), 9 2 (t)), 
that we consider to be the word line (i.e., geodesic) of the object in the coordinate system 
(0i, #i, <p2, #2)- If we see sufficiently many objects, we have sufficiently many geodesies. 

Of course, we cannot get lightlike or spacelike geodesies by this procedure. In the best case, 
we can reconstruct (numerically) sufficiently many geodesies, in the sense their velocity 
vectors are dense in a certain open subset of TM. See also the discussion in [22] . 

Now, as a rule, we can not get the natural parameter (=proper time) of an object. Indeed, 
if the relativistic effects are are not negligible, the proper time of the object is not our own 
time t, i.e., the curve (4>i(t), &i(t), 4>2(t), #2^)) is a reparameterized geodesic only. If we can 
not observe a periodic process on an object (note that the astronomical objects such that 
we can register a periodic process on, for example pulsars, are very rare) or any other way 
to measure the own time of the object, we can not obtain the own time of the objects by 
astronomic observations (see also the discussion in [To]). 

In view of this discussion, the following two problems (Problem [TJ and Problem [2] below) 
in the theory of geodesically equivalent metrics are interesting for general relativity: 

Problem 1. How to reconstruct a metric by its unparameterized geodesies? 
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One can obtain unparameterized geodesies by astronomic observations: 




The general setting is as follows: we have a family of smooth curves j(t] a) in U C M 4 
depending on 6-dimensionafl parameter a = (a±, olq)\ we assume that the family is 
sufficiently big (we formalize 'sufficiently big' in the beginning of Section [2TTj) . We need to 
find a metric g such that for every fixed a the curve j(t; a) is a reparameterized geodesic 
of g. 

Mathematically, the problem has sense in every dimension and for every signature of the 
metric. In dimension 2, versions of this question were considered by S. Lie [M] and R. 
Liouville [3S] , and were also discussed by Veblen and Thomas jlHl SSI E] and Eisenhart 
[13] in the beginning of 20th century. In the realm of general relativity, the problem was 
explicitly stated by J. Ehlers et al [13], where it was said that "We reject clocks as basic 
tools for setting up the space-time geometry and propose ... freely falling particles instead. 
We wish to show how the full space-time geometry can be synthesized ... . Not only the 
measurement of length but also that of time then appears as a derived operation. " 

This problem can be naturally divided in two subproblems. 

Subproblem 1.1. Given a family of curves 7(t;a), how to understand whether these 
curves are reparameterised geodesies of a certain affine connection? How to reconstruct 
this connection effectively? 

We will say that a metric lies in a projective class of a certain symmetric affine connection 
T = Tj k , if every geodesic of g is a reparameterized geodesic of V. 

Subproblem 1.2. Given an affine connection T = T l j k , how to understand whether there 
exists a metric g in the projective class ofT? How to reconstruct this metric effectively? 

Both subproblems were actively discussed in the literature. In dimension 2, the answer on 
Subproblem 11.11 is classical and was known already to Sophus Lie; given a family of curves 

Mocally, the set of unparameterized geodesies of an n— dimensional manifold has the structure of a 
manifold of dimension 2{n — 1) 
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one constructs an ODE of the second order y"(x) = f(x, y(x),y'(x)); the curves ■y(t; a) are 
reparameterized geodesies of a certain connection if and only if the right hand side of the 
ODE is a 3rd degree polynomial in y'(x), 

f(x, y(x), y\x)) = A(x, y(x)) + B(x, y(x))y'(x) + C(x, y(x)) (y'(x)f + D(x, y(x)) (y'(x)f . 

The answer in the multidimensional case can be obtained using the same idea as in dimen- 
sion 2, we give it in Section |2~T1 

The second subproblem is more complicated and is almost open. In dimension 2, the sub- 
problem was considered in the recent paper pjQj of Bryant et al: given an affine connection, 
they construct a system of differential invariants that vanish if and only if there exists a 
metric (in a neighborhood of almost every point) in the projective class of this connection. 
The invariants are very complicated and are of very high orders. 

In theory, one can also obtain a similar answer in every dimension. Indeed, by [12], in 
every dimension the existence of a metric in a projective class is equivalent to the existence 
of a nontrivial solution of a certain overdetermined system of linear PDE in the Cauchy- 
Frobenius form (i.e., the sysem is of first order and all derivatives of unknown functions are 
explicit (linear) expressions in the unknown functions). Given an overdetermined system 
of PDE in the Cauchy-Frobenius form, one can always, in theory, construct a system of 
differential invariants that vanish if and only if the system admits a nontrivial solution 
(in a neighborhood of almost every point). An effective construction of these differential 
invariants could be very complicates. The results of jTU] show that it is indeed the case 
in dimension 2. It is hard to predict whether the system of differential invariants is easier 
in the multidimensional case (normally multidimensional cases are harder than lowdimen- 
sional; but sometimes overdetermined systems are easier to analyse in higher dimensions, 
because they can have higher degree of overdetermination). 

In the present paper, in Section 12.2.21 we give an algorithmic answer to Subproblem 11.21 
under the additional assumption that the metric g we are looking for is Ricci-flat and 
the projective class satisfies certain nondegeneracy assumption, i.e., in a situation most 
interesting from the viewpoint of general relativity. In Section 13.11 we also discuss the case 
of arbitrary metric: we show that also in this case one can algorithmically reconstruct the 
metric by its projective class assuming certain nondegeneracy assumption on the projective 
class; though in this case the nondegeneracy assumption is harder to check. 

Remark 1. Of course it is important in what form the geodesies 7(t; a) are given. Be- 
low, it will be clear what information we need from j(t; a) in order our algorithm works. 
If the geodesies are given numerically (which is the case if they came from astronomic 
observations), this information could be extracted without difficulties. 

Problem 2. In what situations is the reconstruction of a metric by the unparameterised 
geodesies unique (up to the multiplication of the metric by a constant)? 

The example of Lagrange/Beltrami above shows that in certain situations the reconstruc- 
tion is not unique: the geodesies of every metric of constant curvature are straight lines, 
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i.e., the geodesies of the standard flat metric, in a certain coordinate system. Constant 
curvature metrics are not the only metrics that allow nontrivial geodesical equivalence. For 
example, as it was shown by Dini, the following two metrics on U 2 C M? are geodesically 
equivalent 

where X and Y are arbitrary (smooth) functions of the indicated variables such that the 
formulas ([1]) correspond to metrics (i.e., O^X^F ^0 for all (x,y) G U 2 ). This example 
was generalized for all dimensions by Levi-Civita: from his results it is follows that the 
following two 4-dimensional metrics are geodesically equivalent: 

g = (X (x ) - X!(xi))(X (xo) - X 2 (x 2 ))(Xo(xo) - X 3 (x 3 ))dx 2 

+ (X (xo) - X 1 (x 1 ))(X 1 (x 1 ) - X 2 (x 2 ))(X 1 (x 1 ) - X 3 (x 3 ))dx 2 

+ (Xo(x ) - X 2 {x 2 )){X 1 {x 1 ) - X 2 {x 2 )){X 2 {x 2 ) - X 3 (x 3 ))dx 2 1 ) 

+ (Xo(x ) - X 3 (x 3 ))(X 1 (x 1 ) - X 3 (x 3 ))(X 2 (x 2 ) - X 3 (x 3 ))dx 3 3 



9 = XoU) x (xo)x 1 {x 1 )x 2 {x 2 )x 3 (x 3 ) (M*o) ~ X 1 (x 1 ))(X (x ) - X 2 (x 2 ))(X (x ) - X 3 {x 3 ))dx 2 
+ Xl L) Xo(x )x 1 (x 1 )x 2 ( X2 )x 3 (x 3 ) (Mxo) - X 1 {x 1 )){X 1 {x 1 ) - X 2 {x 2 )){X 1 {x 1 ) - X 3 (x 3 ))dx 2 
+ x 2 ( X2 ) Xo(x )x 1 (x 1 )x 2 (x 2 )x 3 (x 3 ) (Xo(x ) - X 2 (x 2 ))(X 1 (x 1 ) - X 2 (x 2 ))(X 2 (x 2 ) - X 3 (x 3 ))dx 2 
+ x 3 U) x (x ) Xl ( Xl )x,(x a )x 3 (x 3 ) (Xo(x ) - XsixzWX^Xi) - X 3 (x 3 ))(X 2 (x 2 ) - X 3 (x 3 ))dx 3 3 . 

(3) 

Here (x , ...,x 3 ) are local coordinates and the functions Xj depend on the indicated variables 
and are such that the metrics have sense. 

In view of this, in the realm of general relativity, Problem |2] can be naturally divided in 
two subproblems. 

We call a metric g geodesically rigid, if every metric g, geodesic equivalent to g, is propor- 
tional to g. 

Subproblem 2.1. What metrics 'interesting' for general relativity are geodesically rigid? 
Subproblem 2.2. Construct all pairs of nonproportional geodesically equivalent metrics. 

Let us comment on these subproblems. The part of the Supbproblem 12.11 that is hard or 
even impossible to formalize is the word "interesting" . Instead of formalizing this notion, 
let us give few results in this direction. 

Probably the metrics that are most interesting in the context of general relativity are 
Ricci-flat nonflat metrics. As it was shown by A. Z. Petrov in [12] (see also [22] and |27j). 

4-dimensional Ricci-flat nonflat metrics of Lorentz signature can not be geodesically equiv- 
alent, unless they are afflnely equivalent 
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(two metrics are affinely equivalent, if their Levi-Civita connections coincide. Affine equiv- 
alent Ricci-flat 4-dimensional metrics are completely understood). It is one of the results 
Petrov obtained in 1972 the Lenin prize, the most important scientific award of the Soviet 
Union, for. 

Recently, the answer of Petrov was generalized in [27] (see also [21]): it was shown that 

if g and g are geodesically equivalent metrics on a 4— dimensional manifold, and g is 
Einstein and of nonconstant curvature, then the metrics are affinely equivalent. 

Let us also give an example of a metric that is important for general relativity and that is 
not geodesically rigid. This is the so-called Friedman-Lemaitre-Robertson- Walker metric 

9 = -dt 2 + R(t) \ d f+^ + 2 d f ; « = +l;0;-l, (4) 
1 + j(x 2 + y 2 + z 2 ) 

where R = R(t) is a real function (the scale factor) of the 'cosmic time' t. The metric is 
not geodesically rigid. Indeed, for every constant c such that the formula below has sense, 
the metric 

-1 , 2 Riff dx 2 + dy 2 + dz 2 
9 (R(t) 2 + c) 2 c( J R(t) 2 + c)l + f(x 2 + y 2 + ^) [) 

is geodesically equivalent to g (one can see it directly as it was done for example [40] or 
[23] . see also discussion in [16]. Actually, the pair of geodesically equivalent metrics (|4|5|) 
is a special case of geodesically equivalent metrics from Levi-Civita [33]). 

For certain functions R, the metric (J4]) is the main ingredient of the so-called Standard 
Model of modern cosmology, and is of cause very interesting for general relativity. 

The metrics listed above, i.e., Einstein metrics and FLRW metrics, are without any doubt 
interesting for general relativity. Of cause, there are other metrics that could be interesting 
for general relativity, and we consider it very important to understand what 'interesting' 
metrics are geodesically rigid. In the present paper, in Section 13. 1\ we prove that almost 
every 4-dimensional metric is geodesically rigid. 

Let us explain what we understand under almost every. Our result is local, so we will work 
in a small neighborhood U C with fixed coordinates (xi, ...,x^ ). We consider a metric g 

n(n+l) n(n+l) 

as the mapping g : U — > M 2 = M 10 ; the space M 2 should be viewed as the space of 
symmetric n x n-matrices. On the space of metrics (viewed as mappings) we consider the 
standard uniform C 2 — topology: the metric g is e— close to the metric g in this topology, 
if the components of g and their first and second derivatives are e— close to that of g. 

In the present paper, we prove that 

for any metric g and every e > there exists a metric g such that g is e-close to g in the 
C 2 — sense, and such that g is geodesically rigid. Moreover, there exists e' > such that 
every metric that is e'— close to g in the C 2 — sense is also geodesically rigid. 

The result is also true in dimensions > 4; the proof is essentially the same. Now, concerning 
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the lower dimensions, the result is true in dimension 3, if we replace the uniform C 2 — 
topology by the uniform C 3 -topology. The proof (will not be given here) is based on the 
same idea. In dimension 2, the result is again true, if we replace the uniform C 2 — topology 
by the uniform C 8 -topology. 

This result was expected, at least if we replace C 2 — topology by C°°-topology. Indeed, by 
Sinjukov [44] and Eastwood et al [32], the existence of a metric geodesically equivalent to 
a given one is equivalent to the existence of a nontrivial solution of a certain linear system 
of partial differential equations in the Cauchy-Frobenius form (1181) . whose coefficients are 
certain invariant expressions in the components of the given metrics and their derivatives. It 
is known that the existence of the solution of such system is equivalent to certain differential 
conditions on coefficients, that is, on the entries of the metrics. If there exists at least one 
metric that is geodesically rigid, then the differential conditions are not identically fulfilled, 
and almost every (in the C°°— sense) metric is geodesically rigid. Now, the existence of 
geodesically rigid metrics in dimensions n > 3 is wellknown (at least since Sinjukov [4"3]). 
The existence of geodesically equivalent metrics in dimension n = 2 is more tricky; it follows 
from Kruglikov [31] where all above mentioned differential conditions were constructed. So 
in a certain sense our result is the improving C°°— closeness (which should be clear to 
experts, though we did not find a place where it is written) to C 2 — closeness. 

Let us now comment on Subproblem 12.21 First of all, the problem is very classical, and 
was explicitly asked by E. Beltrami! in [2]. In the Riemannian case, it was solved by Dini 
in dimension 2 and Levi-Civita in all dimensions. More precisely, Dini has shown that 
locally, in a neighborhood of almost every point of a two-dimensional manifold, every two 
geodesically equivalent metrics are given by the form (ITJ) in a certain coordinate system. 
Levi-Civita has generalized this result to every dimension, we recall his result in Section 

mm 

Unfortunately, the proofs of Dini and Levi-Civita require that the (l,l)-tensor g ggj is 
semi-simple (i.e., has no Jordan blocks), and that all its eigenvalues are real. If one of 
the metrics is Riemannian, this condition is fulfilled automatically. Examples show the 
existence of geodesically equivalent pseudo-Riemannian metrics such that the (l,l)-tensor 
g %i gij is not semisimple or/and its eigenvalues are not real. The examples exist already in 
dimension 2: as it was showrjf] in [6], the metrics from every column of the table 

2 Italian original from [5] : La seconda . . . generalizzazione . . . del nostro problema, vale a dire: riportare 
i punti di una superficie sopra un'altra superficie in modo che alle linee geodetiche della prima corrispondano 
linee geodetiche della seconda 

3 As is was explained in [6] , essential part of the result could be attributed to Darboux [11] 
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Liouville case 


Complex-Liouville case 


Jordan-block case 


9 


(X(x)-Y(y))(dx 2 -dy 2 ) 


^s{h)dxdy 


(1 + xY'(y))dxdy 


9 


( 1 1 \ ( dx 2 dy 2 \ 
\Y(y) X(x) J \X(x) Y(y) J 


( 3W Ydx 2 

\^(h) 2 +K(h) 2 J UX 

1 " M(h)Q{h) dxdv 

i { Q(h) Ydv 2 


1 -^l{-2Y(y)dxdy 
+ (l + xY'(y))dy 2 ) 



are geodesically equivalent (we assume that the functions X and Y depend on the indi- 
cated variables only, and that the function h is a holomorphic function of the complex 
variable z = x + i ■ y). Moreover, every pair of 2-dimensional geodesically equivalent 
pseudo-Riemannian metrics has this form in a neighborhood of almost every point in a 
certain coordinate system. 

By direct calculations we see that the (l,l)-tensor g l£ gej for these metrics is semisimple 
with two real eigenvalues in the Liouville case (we also see that the form of the metrics 
is very similar to ([!]), the only difference is the signature), has two complex-conjugated 
eigenvalues in the Complex-Liouville case, and is not semisimple in the Jordan-block case. 

Actually, certain authors consider that the Subproblem 12.21 is also solved; the solution 
is attributed to Aminova pQ. Unfortunaltely, the author of the present paper does not 
understand her result, and has certain doubts that it is correct. More precisely, in view of 
[TJ Theorem 1.1] and the formulas [U (1.17), (1.18)] for k — 1, n = 4 and all es equal to +1, 
the following two metrics g and g given by the matrices (where u is an arbitrary function 
of the variable x&). 









2X4 












1 




1 




3x 3 + 3u (x 4 ) 2x 2 X\ 



2x 4 - 5 



— X4 



3 £3 + 3 ui (xi) 
2x 2 
■r-i 
4xiX 2 



x 3 +u>(x 4 ) 



—3 23— 3 cj(x4)+2 X2X£ 

3 £3+3 u)(x4)— 2 X2a;4+3;ia;4 2 
1E4 7 



-3x3-3u(x4)+2 3:2X4 3 x 3 +3 u(x 4 )-2 x 2 x 4 +x 1 x 4 2 (-3x s -3uj(x 4 )+2x 2 x 4 )(2x 1 x 4 2 +3x3+3w(x 4 )-2x2X 4 ) 



x 4 ° x 4 ° x 4 ' x 4 a 

should be geodesically equivalent, though they are not (which can be checked by direct 
calculations). Note that the metrics above have signature (2,2), so they are not that in- 
teresting for general relativity. In the case of Lorentz signature, the theorem of Aminova 
seems to be correct, but still it is very complicated to extract the precise formulas from 
her works. 

Note also that according to [T], in the case of Lorentz signature, geodesically equivalent 
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metrics were discribed by Petrov |JT] in dimension 3, by Golikov [18] in dimension 4, and 
by Kruchkovich [30J in all dimensions. From these papers, we were able to find (and to 
check) the paper of Petrov [41] only. 

In the present paper, we combine recent results of {?[/ and above mentioned results of |2J/ 
and to give an easy algorithm how to obtain a list of pairs of all possible geodesically 
equivalent ^-dimensional metrics g,g of Lorentz signature. 

More precisely, we explain (following [7]) that every such pair can be obtained by applying 
the explicit gluing construction from Theorem [3] to building blocks, and provide explicit 
formulas for all possible building blocks. One can easily obtain a complete list of metrics 
by this algorithm. There exists three possible three-dimensional building blocks, three 
possible two-dimensional, and one possible 1-dimensional, so all together there exists 10 
normal forms for geodesically equivalent (nonproportional) metrics of Lorentz signature. 
The normal forms are given by explicit formulas and allow certain freedom as (almost) 
arbitrary choice of functions of one variable or constants or metrics on two- or three- 
dimensional disks. We also explain the (only) difficulty in applying this algorithm in 
higher dimensions. 

2 Problem [D How to reconstruct a metric by its 
unparameterized geodesies. 

2.1 Subproblem [LTD how to reconstruct a connection by unpa- 
rameterized geodesies, and when it is possible. 

We will work in arbitrary dimension n > 2, in a small neighborhood [/ ct" We assume 
that we are given a family of smooth curves ^(t; a). We assume that the family is sufficiently 
big in the sense that at any point x G U the set of vectors 

:= {£ T X0 U | there exists a and t such that ^(l(t; a))^^ is proportional to £} 

contains an open subset of T X0 U. We put Q = [_) xeU £l x - We will call a pair (to; a) 
xo— admissible, if ^ (j(t; a)) , f=to G Q XQ . We need to understand whether there exists a 
symmetric affine connection T such that every curve j(t; a) is a reparameterized geodesic 
of T, and construct this connection if it exists. 

It is well known (at least since the time of Levi-Civita (33]) that, in local coordinates, every 
geodesic 7 : I — > U, 7 : 1 1— > 7 J (t) G U C M n of a symmetric affine connection T is given in 
terms of arbitrary parameter t as solution of 

dt 2 bc dt dt 1 \dt J dt ' 1 ' 
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Better known version of this formula assumes that the parameter is affine (we denote it by 
"s") and reads 

it is easy to check that the change of the parameter s — > t transforms ([7]) in (jH]). 

For further use, let us note that if we linearly change the parameter t of a curve ^(t; a) 
(by putting t = const • t new ), the left hand side of is multiplied by const 2 implying that 
the function / should be homogeneous of degree 1: /(const • £) = const • /(£) for every £ 
(such that £ € f2). This allows us to assume without loss of generality that for every x the 
subset Q x C T X U contains a cone over a nonempty open subset. 

Let us now take a point x G U . For every x — admissible (t ; a), we view the equations (EJ) 
as a system of equations on the entries of r(x ) and on the function /|n ; the coefficients 

in this system come from known data ( ) > [ ~ ^Jt^ I ■ Since we have infinitely 

V / \t = tQ V / \t = tfi 

many xq— admissible (£;a)'s, we have an infinite system of equations. Let us show that 
if this system of equations is solvable, then the solution is unique up to a certain 'gauge' 
freedom. 

Let us first describe the gauge freedom: we consider two connections V and T related by 
Levi-Civita's formula 

ri = Tic - - (8) 

where = 0j is a one form. Suppose the curve 7 satisfies the equation (jHD with a certain 
function /. Substituting Y given by (jSj) in the left hand side of (jSJ) and using 



drdY n (d-f\\d-f a 



dt dt V dt I dt 



we obtain that the same curve 7 satisfies the equation §5§ with respect to the connection 
r and the function 

f(v):=f(v) + 2(v%). (9) 

Thus, if (r, /) is a solution of (j5J), then for every 1— form (f> the pair (F, fj given by (I8"f9~j) 
is also a solution. Let us show that up to this gauge freedom the connection T and the 
function / are unique. 

We again work at one point x £ U and again view (j^D as equations on (r, /). Suppose we 
have two solutions (r, /) and (F, /). We subtract one equation from the other to obtain 

f a bc v b v c = f(v)v a , (10) 

where T = T — T, f = f — f. This equation is fulfilled for all vectors v = v a lying in an 
open nonempty Q X(} C T XQ U. Since the mapping a(u,v) 1— > Yl c u b v c is linear in u and v, it 
satisfies the parallelogram equality 

= a(u + v, u + v) + a(u — v,u — v) — 2a(u, u) — 2a(v, v). (11) 
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Combining (jTTj) with (|T0|) . we obtain 

= f(u + v )(v + u) + f(u - v)(u - v) - 2f{u)u - 2f(v)v 

= (f(u + v) + f(u - v) - 2f{u))u + (/(« + u) - /(« - v) - 2f(v))v. 1 ) 

Taking u and t> to be linearly independent, we obtain 

f(u + v) + f(u-v)-2f(u) = 

f(u + v)-f(u-v)-2f(v) = 1 } 

implying f(u + v) = f{u) + f{y). As we explained above, the functions /, /, and, therefore, 
/, also satisfy const • f{y) = /(const • v ). Then, the restriction of / to a certain nonempty 
open subset Q' C Q xo C T X0 U is linear, i.e., is given by f(v) = 2<f> a v a for a certain 1-form 
<p = (f) a and for all v from Q' . Then, the connection 

ft := f I - <j> b 5 a c - <j> c 8 a h 

has the property that for every (t ; a) such that (^)i t _ to £ ^ the corresponding 7 (t; a) 
satisfies (at t = to) the equation 



dt 2 6c dt dt dt 2 bc dt dt V dt dt 

implying r = T implying that T and T are as in (jSJ) implying / and / are as in ([9]). 

Finally, the connection T and the function /, if they exist, are uniquely determined by the 
unparameterized curves 7 (t; a) up to the gauge freedom 

+ + f^f + 2<t> (14) 

Remark 2. If the function / is linear, i.e., if /(£) = 2(f)^ h for a certain 1— form <p, then, up 
to the gauge freedom, we can take / = 0. Moreover, putting / = we exhaust the gauge 
freedom. 

Let us now explain how to reconstruct the pair (r, /) up to the gauge freedom. We give 
an algorithm how to do it. The algorithm gives also a possibility to understand whether 
there exists such (r, /): we will see it that in order to uniquely reconstruct the (possible) 
entries T(xo)j k of the connection at a point Xq, we will need only finitely many 7 (t; a) 
passing through this point. There exists such (r, /), if for all x the entries of r(x )* fc do 
not depend on the x — admissible (to; a) we used to construct r(x )* fc . 

We will work at a point x ; our goal is to reconstruct the components T(xo) l j k . We take xq- 

admissible (to; a) such that the first component 7^ 0. For this geodesic 7 (to; a), 

we rewrite the equation ([6]) at t = to in the following form: 
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drf_ pi d-y a dy b 
dt ab dt dt 



drf_ pi drf_dr^_ _ dr^_ pn drf_drf_ d 2 -y n drf_ _ drf_ d 2 ^ 1 

dt ab dt dt dt ab dt dt ~ d 2 t dt dt d 2 t ' 

The first equation of (115j) is equivalent to the equation of (EJ) for a = 1 solved with respect 
to / (t?). We obtain the second, third, etc. equations of ( I15j) by substituting the first 
equation of ([TBI) in the equations of (jSJ) corresponding to a = 2, 3, etc. 

We consider now a subsystem of ([TBI) containing the the second, third, etc. equations 
of ([TBI) . We see that the system does not contain the function /. Then, for every Xq- 
admissible (to, a), it is a linear (inhomogeneous) system on the components T(xo)^ k . We 
take a sufficiently big number N and substitute N x ~ admissible generic (to; a)'s in this 
subsystem. 

Remark 3. If n = 4, it is sufficient to take N = 12. We understand the world 'generic 1 
in the following sense: for every n pairs (to, a), the velocity vectors ("^)u_ to are linearly 
independent. 

At every point xo, we obtained an inhomogeneous linear system of equations on n 
unknowns r(x )* fc . 

In the case the solution of this system does not exist (at least at one point xq), there exists 
no connection whose (reparameterized) geodesies are j(t;a). 

If the solution exists at all points, the solution is unique up to the gauge freedom (jSJ). 
Indeed, a solution of the last n — 1 equations of (FIB[) gives us also the values / by the first 
equation of l[TB[) . so the gauge freedom in the equations ( TTB[) is the same as of the equations 
(EJ). Thus, a solution, if it exists, gives us the only up to the gauge freedom candidate for 
the entries F(xo)j k at every point xo such that its geodesies are (reparameterized) curves 
7(t;a). 

Assume now that at every point xq, a solution r(xo))h exists. In order to construct the 
entries T(xo)j k (up to the gauge freedome), we used N xq— admissible curves. In order 
to understand whether all geodesies ^(t; a) are reparameterized geodesies of V, we need to 
substitute all geodesies ^(t; a) in the equation ([6]), and check whether it is fulfilled; in this 
case, it is natural to rewrite the equation in the /—free form 

(d 2 r r gdi b dY\ dr 

V dt 2 bc dt dt J dt 



dj 1 p2 d"/ a d"/ b 
dt ab dt dt 



(d 2 r 

d 2 -y 2 d-y 1 
d 2 t dt 



i pi drf_drf_\ / d-y 1 
<~ 1 ab dt dt 1 1 dt 
_ drf_ dV 
dt d 2 t 



(15) 
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2.2 Subproblem 11.2b given an affine connection T = r* fc , how to 
understand whether there exists a metric g in the projective 
class of T? How to reconstruct this metric effectively? 

2.2.1 General theory. 

We are given a symmetric affine connection T l - k on M n , we need to understand whether 
there exists a metric in the projective class of T. In this section we recall (following 
[HI [12]) the general approach how to do it: the existence of a metric in the projective 
class is equivalent to the existence of a nondegenerate solution of a certain system of linear 
PDE in the Cauchy-Frobenius form, and, in theory, there exists an algorithmic way to 
understand the existence of such solutions. 

Theorem 1 ([12], see also references inside), g lies in a projective class of a connection 
T l j k if and only if a ab := g ab ■ det(g) 1//(n+1 ) is a solution of 



(V^) - ^ (V^S: + V t aS b ) = 0. (16) 

Here a ab := g ab ■ det(#) 1/(n+1) should be understood as an element of S 2 M <g> (K n ) 2 ^ n+l ^M. 

In particular, V a a bc = ^a bc + T b ad a dc + T da a bd - -^—T d da a bc 

ox a n + 1 

" v — ' » w ' 

Usual covariant derivative addition coming from volume form 



The equations ( Fl6|) is a system of ^!Li^±l) _ n \ linear PDEs of the first order on n ^" 2 +1 ^ 
unknown components of a. 

Two-dimensional version of these equations was essentially known to R. Liouville [35]: 
instead of working with a ab := g ab ■ det(g) 1 ^ n+1 \ he worked with a^- = det A 2/3 ^j; in 

dimension 2 the entries of <r u and are linearly related. The 2-dimensional analog of the 
equations (fl6|) is then the Liouville system of 4 PDE's of the first order 



(17) 



+ 2 K a 12 - 2/3 K ian = 

+ § W + 2 K oa 2 2 + 2/3 K x a X2 - 4/3 K 2 a u = 

^ + 2^ + 4/3^022-2/3^2- 2 K 3 au = 

^ + 2/3 K 2 a 22 - 2 K 3 a 12 = 0, 

where K := -T 2 n , K x := - 2r? 2 , K 2 := -F 2 22 + 2T\ 2 , K 3 := Y\ 2 . 

Remark 4. One sees that the gauge freedom ([Ml) does not affect the coefficients Kq, K 3 
of the equation (IT7|) . One can check by calculations that this is also true in all dimensions: 
the gauge freedom (fT4"]) does not change the equations ffTBT) . 

The PDE-system (fl6|) can be prolonged (see [l2j) to the system 
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V a <x bc 



<*a V + S a C fM b 



cd 



(18) 



where P is the symmeterized Ricci-tensor, y the Cotton- York- Tensor and W a b c d the pro- 
jective Weyl tensor for the connection T. 

Remark 5. Here we use another index convention for the projective Weyl tensor than in 
Section 12.2.21 of our paper. This convention is the same as in [12] , and is standard in the 
so-called tractor calculus, we refer to [12] for precise formulas. In Section 12.2.21 we will 
explain the convention used there by given the formula for Weyl tensor. 

The system (fTSj) is a linear system of PDE of the first order on the unknown functions a bc , fi b 
, p. Moreover, all derivatives of unknowns are expressed as functions of unknowns, i.e., the 
system is in the Cauchy-Frobenius form. One can understood this system geometrically as a 
connection on the projective tractor bundle = £ l > bc \-2)+£ b (-2) + £(-2), see [12] for 

details. The solutions of the system are then parallel sections of the connection; there exists 
an algorithmic way to understand whether a certain connection admits a nontrivial parallel 
section. In the two-dimensional case, the algorithm was fulfilled for certain projectively 
homogeneous connections in [9]; for arbitrary two-dimensional connection, the algorithm 
was fulfilled in [10], and the answer (i.e., the differential conditions on K~i such that its 
vanishing implies the existence of a nontrivial solution) appears to be very complicated. 
In theory, one can fulfill this algorithm for every dimension; it is clearly a nontrivial task. 
In the next section we will show that, under the additional assumption that the searched 
metric is Ricci-flat, there exists a trick that simplifies the algorithm. 

2.2.2 The case n = 4, g is Ricci-flat. 

Let us now assume that we know the geodesies of a nonflat Ricci-flat metric. That is, we 
know a certain T such that for a certain <p a which we do not know T% c := T% c + 5%4> c + ^c^c 
is the Levi-Civita connection of a certain nonflat Ricci-flat metric which we again do not 
know. Our goal is to find this metric (which I call g). By the above mentioned results of 
Petrov [JT] j Hall et al [221 121] , and Kiosak et al [27] , the metric is unique up to multiplication 
by a constant; the goal of this section is to explain how to find it algorithmically. The 
algorithm works under certain additional (generic) condition on the connection T. 

We consider the projective Weyl tensor introduced in [38] (not to be confused with the 
conformal Weyl tensor) 



(in our convention Rj k = R a jk a , so that W a jk a = 0). 

Weyl has shown that the projective Weyl tensor does not depend of the choice of connection 
within the projective class: if the connections T and T are related by the formula (jSJ), then 



W l jkt :- R l jk £ - — 



(5\ R jk - SI Rji) 



(19) 
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their projective Weyl tensors coincide. Now, from the formula (1191) . we know that, if the 
searched g is Ricci-flat, projective Weyl tensor coincides with the Riemann tensor R l j ke of 
g. Thus, if we know the projective class of the Ricci-flat metric g, we know its Riemann 
tensor. 

Then, the metric g must satisfy the following system of equations due to the symmetries 
of the Riemann tensor: 

f giaW a jkm + gjaW a ikm = 

1 g ia W a jkm - g ka W a mij = 1 > 

The first portion of the eqautions (120]) is due to the symmetry (R^km = —Rjikm), and the 
second portion is due to the symmetry (Rkmij = Rijkm) of the curvature tensor of g. 

We see that for every point Xq E U (120]) is a system of linear equations on g(xo)ij. The 
number of equations (around 100) is much bigger than the number of unknowns (which 
is 10). It is expected therefore, that a generic projective Weyl tensor W l ju admits no 
more than one- dimensional space of solutions (by assumtions, our W admits at least one- 
dimensional space of solutions). The expectation is true, as the following classical result 
shows 

Theorem 2 (@2l US 120] [ZD [38]). Let W* : jke be a tensor in M 4 such that it is skew- 
symmetric with respect to k, I and such that its traces W a a M and W a j a e vanish. Assume 
that for all 1-forms ^ we have W a jkeC,a ^ 0. Then, the equations ( 120]) have no more 
than one- dimensional space of solutions. 

Let us comment on the condition W a jke^, a ^ 0. In this context, for every fixed indexes 
k, £, W l j^ could be viewed as a n x n-matrix; and the condition W a j^ a = means that 
the matrix has a nontrivial kernel (in particular, it is degenerate). Now, the condition 
W a jke£, a = means that for all indexes k, I the kerns of the n x n-matrices W a jke have 
nontrivial intersection. Thus, it is a very restrictive condition on W , and, therefore, on 

r. 

This result shows that, under the assumptions that for all ^ =^ we have W a jke^, a ^ 0, 
we can reconstruct the conformal class of the metric g by solving the system of linear 
equations (12"U]) . This can be done algorithmically. Then, we also know the conformal class 
of a in f ]TB]) . i.e., we know that a is of the form 

a lj = eV, (21) 

where a u is known and comes from the solution of the linear system (120]) . and the function 
A is unknown. Substituting the ansatz (T2T]) in the system f]T6]) . we obtain an inhomogeneous 
system of linear equations on the components J^. Direct calculations show that this system 
has at most one solution; since we assumed the existence of the metric in the projective 
class, one can always solve this system and obtain all Finally, we can obtain the 
function A, and, therefore, the metric g, by integration. 
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Let us note that in all steps we assumed that a Ricci-flat metric g exists in the given 
projective class. But the algorithm also gives us an algorithmic check whether such metric 
exists: one should go along the steps of the algorithm and look whether something goes 
wrong. 

For example, the system ( 120]) could have no nontrivial solution (i.e., every solution g^ of 
(j20l) has zero determinant). Then, no Ricci-flat metric g exists in our projective class. 

If the system ([2D]) has nontrivial solution, then, after plugging the ansatz (f2~T]) in (fl6l) . we 
obtain a system of nonhomogeneous linear equations on This system may have no 
solution at all (the number of equations is much bigger than the number of unknowns; 
besides, the system is inhomogeneous) , or the 1— form -§^dxi may be not closed. In this 
case, no Ricci-flat metric g exists in our projective class. 

Finally, if the system (j20l) has nontrivial solution, if we can solve the system of linear 
equaitons we obtain after plugging the ansatz ( f2TT) in ( ITB"]) . and if the solution satisfies 
the 'closeness' condition = tttttt, then we do obtain a metric qa in the projective 

ax K ox 1 ox 1 ax K ' f'J j 

class. The metric must not be Ricci-flat though. 

Remark 6. In Section I3TT] we show that one can reconstruct an almost every (4-dimensional) 
metric by its projective class, see Remark [9] there. In the case of arbitrary metric, the non- 
degeneracy assumption on the projective class is more complicated, and it is harder to 
check it. 

3 Problem 2: In what situations is the reconstruc- 
tion of the metric by the unparameterised geodesies 
unique (up to the multiplication of the metric by a 
constant)? 

3.1 For generic 4-dimensional metric, the reconstruction of the 
metric by the unparameterized geodesies is unique. 

Let us first construct one geodesically rigid metric in dimension n — 4. 

Using the formula (jT9j) . by short tensor calculations we see that the metric g^ must satisfy 
the equation 

ng a{i Wi\ kl = g ab W {i ab[l 5% (22) 

where n = 4, the brackets "[ ]" denote the skew-symmetrization without division, and the 
brackets "( )" denote the symmetrization without division. 

Remark 7. Actually, the equation f[2"2"]) is a part of the curvature of the tractor connection 
( [18]) ; in this context, it was obtained in [T2j . 
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V 



We take a 4-dimensional metric g such that at the point x$ it is given by the identity 
matrix 

/I \ 

1 

1 

V 

and such that its curvature tensor (with lowered indexes) Rijki at the point Xq is given by 

Rijki = hikhji — huhjk + H ik Hji — HuHj k , (23) 
where the entries at xq of the (0, 2)— tensors h and H are given by the diagonal matrices 





h 



1 







, H 







1 



Such metric g exists by p21 Theorem 1.12.2] (see also [HI Theorem 1.1]), since the tensor 
( 1231) satisfies all symmetries of the curvature tensor. 

Every metric g geodesically equivalent to g has the same projective Weyl tensor as g. We 
view the equation ( 1221 as the system of homogeneous linear equations on the components 
of g\ every metric g geodesically equivalent to g satisfies this system of equations (with the 
same coefficients W\). At the point xo, this is a system on 10 unknowns g(xoY J . Since the 
system is symmetric in i,j and skew-symmetric in k,l, the system contains 60 equations 
(actually, less because of certain hidden symmetries inside). By direct calculations, we see 
that the rank of this system is 9. Indeed, it has at least one nontrivial solution, namely 
g(xQ)^ , so its rank is at most 9. One can easily find 9 linear independent equations of this 
system (so the rank is at least 9), namely the equations corresponding to the followings 
indexes k, I): 



{i,j,k,l) 


equation 


(1,1,2,1) 


-12# 12 = 


(1,1,3,1) 


2g 13 = 


(1,1,4,1) 


2g u = 


(2,1,2,1) 


5g n - Qg 22 + g 33 = 


(2,1,4,1) 


g 24 = o 


(2,2,3,2) 


8g 23 = 


(3,1,3,1) 


-Ag 11 + g 33 + Ag 22 - g u = 


(3,1,4,1) 


2g 34 = 


(3,2,3,2) 


-Qg 22 + Ag 33 + 3g u - g u = 0. 



We see that the equations in the table are linearly independent. Thus, at the point Xo, the 
set of solutions of this system is 1-dimensional, implying that every metric g, geodesically 
equivalent to g, is proportional to g. 
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Let us show that at every point in a small neighborhood of xo, the system (122]) also has 
rank 9. Indeed, the rank of a matrix is the biggest dimension of a nondegenerate quadratic 
submatrix and therefore is a lower semi-continuous (integer valued) function, i.e., rank 
of this system is at least 9 at every point of a small neighborhood of xq. Now, at every 
point the components g^ give us a nontrivial solution, so the rank can not be bigger than 
9. Thus, in a small neighborhood of Xq, every metric g geodesically equivalent to g is 
conformally equivalent to g. Now, by Weyl |48j, two conformally equivalent 4-dimensional 
metrics are proportional. Then, the metric g is geodesically rigid. 

Now let g be an arbitrary metric in a small neighborhood of xq. We consider the metric 

9t ■■= (1 - t)g + tg. 

The system ( 122]) constructed for this metric has rank 9 for t lying in a small interval around 

1. Since the coefficients of the system are algebraic expressions in t whose coefficients are 
algebraic expressions in the components of g, g and their first and second derivatives, for 
almost all t the system (122]) constructed for the metric g t has rank 9. We take t close to 
such that the metric gt is e— close to g and such that the system ( 122]) constructed for 
the metric g t has rank 9. As we explained above, this metric is geodesically rigid. Every 
metric g that is C 2 — close to gt is also geodesically rigid, since the entries of W for g are 
algebraic expressions in the components of g and its first and second derivatives. Hence, 
the coefficients in the system (122]) constructed for g are close to that of the system (122]) 
constructed for g t implying the system also has rank 9 implying the metric g is geodesically 
rigid as well. 

Thus, for every 4- dimensional metric g and for any e > there exists a metric g t that is 
e— close in the C 2 — sense to g and e' > such that all metrics e' — close in the C 2 — sense 
to gt are geodesically rigid. 

Remark 8. As we mentioned in the introduction, a similar proof can be done for all di- 
mensions n > 4. For dimensions 2 and 3, the proof does not work anymore, since the 
system (122]) has corank at least 2 for all metrics g (one can prove it using the methods of 
[28] §2.3.2]). One can still modify the proof replacing the system (122]) by another projec- 
tively invariant system of equations. This other projectively invariant system of equations 
requires higher derivatives of the components of g though. In dimension 3, one can con- 
struct (using the curvature of the tractor connection (Tl8]) . see also [39]) such projectively 
invariant system such that its coefficients depend on the components of the metrics and its 
first, second and third derivatives. Therefore, for every 3-dimensional local metric g and 
for any e > there exists a metric g t that is e— close in the C 3 — sense to g and e' > such 
that all metrics e'— close in the C 3 — sense to g t are geodesically rigid. Now, in dimension 

2, the construction of the projectively invariant system is much more involving (see |lUj ) 
and requires 8 derivatives of the components of the metric. 

Remark 9. We also see that the projective class of almost every (in the C 2 — sense) 4- 
dimensional metric determines its conformal class uniquely: one can find the conformal 
class by solving the system (122]) . Then, one can proceed along the algorithm from Section 
12.2.21 and understand whether there exists a metric in the projective class, and find it. 
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3.2 Normal forms for pairs of geodesically equivalent 4-dimensional 
metrics such that one of them has Lorentz signature. 



3.2.1 Splitting and gluing constructions from [7]. 

Given two metrics g and g on the same manifold, we consider the (1, 1)— tensor L = L(g,g) 
defined by 

where g lk is the contravariant inverse of g^. 

Remark 10. If n is even, the tensor L is always well defined. If n is odd, the ratio 
det(g)/det(g) may be negative, and the formula (12"4"|) may have no sense. In this case, 
we replace g by — g and make the ratio det(g)/det(g) positive and L well defined. In the 
cases interesting in our context, g and g have the same signature, and the problem with 
the sign does not appear at all. 

Remark 11. The tensor L l - defined in (j24p is essentially the same as as the tensor introduced 
by Sinjukov (see equations (32, 34) on the page 134 of the book jH], and also Theorem 
4 on page 135) and which is often denoted by tensor in the related literature. More 
precisely, L* = aejg h . It is also closely related to a from §2.2.11 g is geodesically equivalent 
to g, if and only if a ab := L a e g eb ■ det(#) 1/(n+1) is a solution of (flB]) . 

The simplified version of the gluing construction does the following. Consider two manifolds 
Mi and M 2 with pairs of geodesically equivalent metrics hi ~ hi on Mi and hi ~ h 2 on 
M 2 . Assume that the corresponding (1, l)-tensor fields Li = L(hi, hi) and L 2 = L(h 2 , h 2 ) 
have no common eigenvalues in the sense that for any two points x G Mi, y G M 2 we have 

Spectrum Li(x) D Spectrum L 2 {y) = 0. 

Then one can naturally construct a pair of geodesically equivalent metrics g ~ on the 
direct product M = Mi x M 2 . These new metrics g and (7 differ from the direct product 
metrics hi + h 2 and hi + h 2 on Mi x M 2 and are given by the following formulas involving Li 
and L 2 : we denote by Xii * = 1, 2, the characteristic polynomial of L^. x% — det(t • 1 — Li). 
We treat the (1, 1)— tensors Lj as linear operators acting on TMj. A polynomial /(L) in L 
is then the (1, l)-tensor of the form f(L) = a (x) ■ Id + ai(x)L + a 2 (x)L 2 + ■ — h a m (x)L m . 
For two tangent vectors 



we put 



u — ( ui , u 2 ) , v — ( Vi , v 2 ) G TM 
gtMi stm 2 eTMi eTM 2 

^(u,u) = hi(x2(L 1 )(ui),vi) + h 2 (xi(L 2 )(u 2 ),v 2 ) , (25) 

= —7Xshi(X2(Li)(ui),Vi) + — y-rhi (xi(L 2 )(u 2 ), v 2 ) . (26) 
X2{0) Xi(0) 
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The corresponding (1, 1)— tensor L = L(g,g) is the direct sum of L\ and L 2 in the natural 
sense: for every 



It might be convenient to understand the formulas (1231 [261) in matrix notation: we consider 
the coordinate system (a; 1 , x r , y r+1 , y n ) on M such that a;— coordinates are coordinates 
on M\ and y— coordinates are coordinates on M 2 . Then, in this coordinate system, the 
matrices of g and g have the block diagonal form 



Theorem 3 (Gluing Lemma from [7]). If hi is geodesically equivalent to hi, and h 2 is 

geodesically equivalent to h 2 , then the metrics g,g given by ( 1231 1231) are geodesically equiv- 
alent too. 

The splitting construction is the inverse operation. We will not describe it completely (and 
refer to [7]); we will use its following corollary explained in [7J §2.1]: 

Every pair of geodesically equivalent metrics h and h in a neighborhood of almost every 
point can be obtained (up to a coordinate change) by applying splitting construction to 
building blocks. 

By a building block we understand an open neighborhood U C lR m with a pair of geodesically 
equivalent metrics h ~ h such that at every point the tensor L given by (12^1) has only 
one real eigenvalue, or two complex-conjugate eigenvalues, and such that the geometric 
multiplicity of the eigenvalue is constant on U. 

Remark 12. Riemannian version of the splitting/gluing constructions was known before, 
see for example [371 Lemma 2] and [361 §§2-2, 2.3]. 

Example 1 . In the definition of the building block, we allow the dimension m = 1. Then, the 
following two metrics on the interval Id 1 with the following two geodesically equivalent 
metrics h = dx 2 and h = X(x)dx 2 (where the function X never vanishes) form a building 
block. Actually, up to a coordinate change, (Ui,h,h) is the only 1-dimensional building 
block. 

Example 2. All possible examples of two-dimensional building blocks can be extracted 
from the table of 2-dimensional geodesically equivalent metrics from the introduction. The 
metrics from the first column of the table do not correspond to a building block, since the 
tensor L for these metrics has two different eigenvalues, X(x) and Y(y). But the metrics 
from the second and the third columns do correspond to the building block, since the 
tensors L for these metrics are given by the matrices 



v = ( vi , v 2 )e T {Xiy) {Mi x M 2 ) we have L{£) = (L 1 (v 1 ),L 2 (v 2 )) • 



£T x Mi &T y M 2 




(27) 



%h) 



Y (x 2 ) 
l + Xl £Y(x 2 ) Y(x 2 ) 
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Of cause, in every dimension, in particular in dimension two, there exists a trivial building 
block (U, h,h = const • h); the tensor L for this metric is a multiple of <5*. From the results 
of [6] it follows that every two-dimensional building block has one of these three forms. 



The formulas for the 3-dimensional building block can be obtained using Petrov [31] and 
Eisenhart [15]; we will give them later. From linear algebra it follows that if the metrics 
g, g have Lorentz signature, then 4-dimensional building blocks are not possible (except 
for the trivial block corresponding to proportional metrics g ~ g := const • g), since in the 
Lorentz signature a g-selfadjoint (1, l)tensor L can not have a Jordan block of dimension 
> 4 with real eigenvalue, and a Jordan block of dimension > 2 with complex eigenvalue. 

Example 3 (Dini formulas (pjj follow from splitting-gluing constructions.). We consider the 
two 1-dimensional building blocks 

h,h x = dx 2 , hi = x ^ x y dx2 ^j and ( J 2, h 2 = -dy 2 , h 2 = ~ y^ 2 d 2/ 2 

We assume that X(x) > Y(y) for all (x, y). The corresponding tensors L\ and L 2 (we view 
them as 1 x 1-matrices) and their characteristic polynomials are 

Li = (X(a;)); L 2 = (Y{y))- Xl (t) = t - X(x) ; X2 (t) = t - Y(y). 

We see that the metrics hi, h 2 satisfy the assumptions in Theorem [3j Plugging these 
data in the formulas ( 1271) . we obtain geodesically equivalent metrics g and g given by the 
matrices 

X{x)-Y{y) \ _ _ ( X {x)*Y(y) 



X{x) - Y(y)) • y ~ y |M=^. 

We see that these metrics are precisely the Dini metrics (pj). For further use let us note 

/X(x) 

that the tensor fT24l) for these metrics is given by L 



Y(y) / 

Example 4 (Levi-Civita metrics ( I2"f3"|) follow from splitting-gluing constructions.). We take 
4 pairs of geodesically equivalent metrics on the interval /. 

9i = dx\ ~ gi = Xl ( Xl) a dxl ; g 2 = -dx\ ~ g 2 = - x ^ x2? dx\ ; 

^3 = dx\ ~ ^ 3 = X3 ( X3)2 dxl ; 9a = ~dx\ ~g A = - x ^y dx\. (28) 
We assume that for i ^ j Xj(xj) ^ Xj(xj) for all x iy Xj G /. 

Gluing (I,g 1 ,g 1 ) and (I,g 2 ,g 2 ), ((I,g 3 ,g 3 ) and (I, g^ g 4 ), respectively) we obtain two pairs 
of geodesically equivalent metrics (we denote them by hi ~ hi {h 2 ~ h 2 , respectively)) on 
the two-dimensional disk U 2 — I x I. These metrics and the corresponding tensors ([21 
were essentially constructed in Example [3] and are given by matrices 

, (Xl(xi) - X 2 (X 2 ) \ -r _ ( Xi(xi)*X2(x2) 
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We see that the metrics hi, ti2 satisfy the assumptions in Theorem [3j Gluing these metrics, 
we obtain the metrics (12|3|) . 

Remark 13. By changing the sign of the metrics fl4]) we can make geodesically equivalent 
metrics g ~ g of arbitrary signature. 

Example 5 (General Levi-Civita metrics). We take m building blocks: the first r building 
blocks are 1-dimensional, and the last m — r building blocks h r+ \ ~ h r +i, h m ~ h m have 
dimensions fcj > 2, i — r + 1, ...,m — r. For cosmetic reasons we think that the first r 
building blocks are 

(Ul, hi = ±dxi 2 , h~i = ± —- . 2 dxi 2 ) , i = l,...,r, (29) 

A-i \%i ) 

the sign ± in hi and /ij is the same for each i, but may be different for different z's. The 
last m — r building blocks are 



U- l ,hi= jr (h i (x i )) ai f) i dx? i dxf i ,hi = -^7 (h^x^a^dx^dx^ 



i = r + 1, m. 



Here the functions Xj are constant for i > r and depend only on the corresponding variable 
Xi for % <r. As above, we assume that Image(Xj) n Image (Xj) = for % 7^ j. The metrics 
hi, i = r+1, m can be arbitrary, but their entries {hi) a .p. must depend on the coordinates 



Xi = (xj, ...,Xt l ) only. 



Inductively applying the gluing procedure, we obtain for g and g the following form: 

(30) 



where 

fl==±n<*-*A « : =^rfx- (31) 

(the signs ± in ( 13T]) depend on the choice of the signs ± in ( 129]) and can be arbitrary). This 
is precisely Levi-Civita's normal form for geodesically equivalent (Riemannian) metrics 
from [331. 



Now, since every pair of geodesically equivalent metrics (in a neighborhood of almost every 
point) can be obtained by a gluing construction, and since in the Riemannian signature only 
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the blocks used above can be used, every Riemannian geodesically equivalent metrics have 
the form (|30|) in a certain coordinate system. This is the famous Levi-Civita's Theorem 
from [33] . 

Note, than the Lorentz signature of g and g does not allow the tensor L to have complex 
eigenvalues of algebraic multiplicity greater than one. Similarly, it does not allow the tensor 
L to have a Jordan block of dimension 4, or two Jordan blocks. Thus, in order to obtain 
the description of nonpropotional 4-dimensional geodesically equivalent metrics of Lorentz 
signature, one needs the building blocks of dimensions 1,2,3 only. In dimension 1, only 
one buiding block, namely the one from Example [U is possible. 

Geodesically equivalent metrics such that the tensor L has the 2-dimensional Jordan-block 
structure 

'A 1\ (X \ ( a 
XJ ' V X) ' \-(3 a 

were described in Example [2j For the Jordan-block structure 



(32) 




the description of the metrics follows from Petrov [41]: the metrics are given by 
9 = ( 4 x 2 (^A (z 3 )) + 2 ) dx x dx 3 + dx 2 2 



9 



2 (a a ( x3 dx2dx3 + xi * (dt A ( x3 ^ 

4 x 2 A (x 3 ) 2 A (x 3 )) + 2 A (x 3 ) 2 ) d Xl dx 3 + A {x 3 ] 



2 cfe 2 2 



(4x 2 A (x 3 ) (^A (x 3 )) + 2A (x 3 ) - 2x x A (x 3 ) 2 (^A (x 3 ))) dx 2 dx 3 
' 4 ^ {lk X M) 2 + 4 ^ (afc A (x 3 )) - 4 W (x 3 ) A (x 3 )) cfe 3 



(33) 



+ [l+ Xl *\( X3 f(jL\( X3 )) -2x 1 X{x 3 )[£X(x 3 )))dx 



d_ 

d,X3 



The corresponding L is given by the matrix 



(afeM**))*! 



A (x 3 ) 1 

A(x 3 ) 2 (-^-A (x 3 )) x 2 + 1 
A(x 3 ) 











Remark 14. Actually, the formulas (133]) are slightly more general than that of j4T]. They 
are equivalent to the formulas from j4T] (modulo a coordinate transformation) at the points 
such that dX ^ 0. The formulas [4T] were obtained together with A. Bolsinov; they can be 
generalized for every dimension. We will publish this result elsewhere. 
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As it follows from J2?l Lemma 6], if L has the Jordan-form 



1 

A 



the eigenvalue 



A, 



A is constant, and the metrics are affinely equivalent (i.e., Levi-Civita connections of g 
and g coincide). Affinely equivalent metrics whose tensor L has this form were essentially 
described by Eisenhart in [T5], see also [2HJ Theorem 1]. From their description it follows, 
that, in a certain coordinate system, geodesically equivalent metrics g ~ g are given by 

g = 2 dx 3 dxi + h (x 2 , x 3 ) n dx 2 2 + 2 h(x 2 , x 3 ) 12 dx 2 dx 3 + h(x 2 , x 3 ) 22 dx 3 2 , 
g = 2 a dx 3 dxi + ah (x 2 , x 3 ) u dx 2 2 + 2 ah (x 2 , x 3 ) 12 dx 2 dx 3 + /3dx 3 2 + ah (x 2 , x 3 ) 22 dx 3 2 , 

(34) 



where a and (3 are constants. 
Now, the metrics g,g such that L 



'A 



A are conformally equivalent. By by the 
\ 

classical result of Weyl [H], they are proportional (i.e., g = const • g). 

Thus, we have described all building blocks that can be used in constructing metrics of 
Lorentz signature; Theorem [3] gives us the construction. Let us count the number of cases 
in dimension 4: we can represent 4 as the sum of natural numbers by 4 different ways: 



Dim of blocks 


Description of blocks 


# of cases 


1+1+1+1 


All building blocks are as in Example |1J and g ~ g are 
essentially fT2|3]) with the changed sign of dx\ 


1 


1 + 1+2 


The first two building blocks are as in Example m the 
third is as in Example |2] 


3 


2+2 


Both building blocks are as in Example |2J at least one 


3 


1+3 


Tlie first building block as is Example [1| the second is 
as in (|33j). as in (|34l). or trivial 


3 



Remark 15. The general schema also works in higher dimensions, but in this case there is 
the following essential difficulty (and this is the only difficulty): up to our knowledge, for 
dimensions n — 1 > 5, there is no description of all pairs of (g, L) such that g has Lorentz 
signature and L is an (l,l)-selfadjoint tensor such that it is covariantly constant, and such 
that the Jordan normal form of L is 



/A 



1 

A 1 
A 



\ 



V 



A 
V 



(35) 



In dimension n = 4, since n — 1 = 3, the Jordan normal form (13"o*|) coincides with (I32p . and 
the local description follows from [H]. In dimension n = 5 we have n — 1 = 4 and one can 
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obtain the local description (we will not do it in the present paper) combining the results 
of [151 122] with the algebraic description of possible holonomy groups of 4- dimensional 
metrics of Lorentz signature (see e.g. [25l EE])- 
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